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Abstract. The explicit semiclassical treatment of logarithmic perturbation theory 
for the nonrelativistic bound states problem is developed. Based upon ?i-expansions 
and suitable quantization conditions a new procedure for deriving perturbation 
expansions for the one-dimensional anharmonic oscillator is offered. Avoiding 
disadvantages of the standard approach, new handy recursion formulae with the same 
simple form both for ground and exited states have been obtained. As an example, the 
perturbation expansions for the energy eigenvalues of the harmonic oscillator perturbed 
by Xx^ are considered. 



1. Introduction 

Logarithmic perturbation theory [1-8] is one of the principal approximation techniques 
in theoretical and mathematical physics. Within the framework of this theory, the 
conventional way to solve a quantum-mechanical bound-state problem consists in 
changing from the wave function to its logarithmic derivative and converting the time- 
independent Schrodinger equation into the nonlinear Riccati equation. Such a procedure 
leads to handy recursion relations in the case of ground states, but becomes extremely 
cumbersome in the description of radial excitations when nodes of wave function are 
taken into account. Although several attempts have been made to improve the method 
in the latter case [9-11], they have not resulted in a desirable simple algorithm. 

On the other hand, it is well known, that the radial quantum number, n, 
most conveniently and naturally is introduced into consideration by means of 
quantization conditions, as in the Wentzel-Kramers-Brillouin (WKB) approximation 
[12-14]. However, the WKB approach is more suitable for obtaining energy eigenvalues in 
the limiting case of large quantum numbers, whereas the perturbation theory deals with 
low-lying levels. Usually, the perturbation results are obtained within the framework of 
the WKB method by recasting the WKB expansions [15-19]. 

The objective of this paper is to develop an explicit semiclassical treatment 
of logarithmic perturbation theory and to describe a straightforward semiclassical 
procedure for obtaining the perturbation corrections through handy recursion formulae, 
having the same form both for ground and excited states. 

For the sake of simplicity, we restrict ourselves to the consideration of the bound- 
state problem for the one- dimensional anharmonic oscillator. The generalization to 
the three-dimensional problem reqiures including the quantization conditions for orbital 
momentum and will be published elsewhere. 



2. Method 

The system to be treated is described by the Schrodinger equation 

-—U"{x) + V{x)U{x) = EU{x) (1) 

where the potential function, V{x) , has a simple minimum and hence can be given by 
the expression 

V{x) = \muj^x^ + Y,hx^+\ (2) 

With changing the scale of the variable, x -^ \/hx , it becomes obvious that 
the coupling constants, /« , appear in common with powers of Planck's constant, h 
. Therefore, the perturbation series must be semiclassical ^-expansions, too. Although 
well known in the folk wisdom of theoretical physics, this assertion was, nevertheless, 
proved not so long ago [20]. It has been argued that the energy eigenvalues under 
consideration should be concentrated near the minimum of the potential and should 
behaved as 

E = nuj{n + \)+Y,Ei{u,n)n\ (3) 

To our knowledge, there are only a few procedures for computing the coefficients 
Ei{uj,n). They involve, in particular, applying the methods of the comparison equation 
[21] and complex "sprout" [22]; an analytic continuation in the h -plane [23]; various 
approaches within the framework of the WKB approximation [15-19]; quantization using 
the methods of classical mechanics [17,24]; and, lastly, expansions in the h ' -series [25]. 
However, all of these methods have some disadvantages. 

Here we propose a new, simpler and more straightforward semiclassical technique. 
Being explicitly opposed to the WKB approach, it is based on different quantization 
conditions which are more appropriate for describing the solution of the bound-state 
problem in the vicinity of a potential minimum. 

Following usual practice, we apply the substitution, C{x) = hU'{x)/U{x) , accepted 
in the logarithmic perturbation theory and go over from the Schrodinger equation (1) 
to the Riccati equation 

hC'{x) + C\x)=2m[V{x)-E]. (4) 

We attempt to solve it in a semiclassical manner with series expansions in the Planck 
constant 

oo oo 

E = J2 Ekh' C{x) = J2 Ck{x)n^ (5) 

k=Q k=Q 
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that result in the system 

Cl = 2m[Vix) - Eo] 

C'q + 2Co{x)Ci{x) = -2mEi 

(6) 

k 

Ci_,{x)+J2a{x)Ck-^{x) = -2mEk. 



i=0 



In the case of ground states, this system coincides with one derived by means of the 
standard technique and can be solved straightforwardly. However, complications of the 
logarithmic perturbation theory arise in the description of radial excitations when the 
nodes of wavefunctions are included in some separate factor. We intend to circumvent 
these difficulties by making use of the quantization conditions. The matter of the latter 
consists in applying the principle of argument, known from the analysis of complex 
variables, to the logarithmic derivative, C{x) . 

Since the wavefunction of the nth radially excited state has n real zeros we have 

1 /■ 1 °° /■ 

— j Cix) dx = —J2^'"f Ck{x) dx = hn. (7) 

fc=0 

There is, however, one important point to note. In the WKB approach, this condition 
is supplemented by the following rule of achieving a classical limit: 

h ^ n— ^oo hn = const (8) 

accompanied by the equality of the quantum energy to the classical one. 

In contrast to this, our method, dealing with low-lying states and being 
complementary to the WKB approach, involves the alternative possibility 

h ^ n = const hn ^ (9) 

which was formerly applied in deriving coefficients of the 1/A^ expansions [26-29]. In 
the limiting case, as ^ — *> 0, a particle is now lowered to the bottom of a potential well 
and its classical energy becomes Ed = minV{x) , which equals to zero in our case. 
Thus, in view of the rule (9), the quantization conditions (7) become 

(fCi(x)dx = n — (fCk(x) dx = k > 1. (10) 

27ri J 2m J 

A further application of the theorem of residues to the explicit form of functions 
Ck{x) easily solves the problem of taking into account nodes of the wavefunctions. 



3. Recursion formulae 

Let us consider tlie system (6) and investigate the behaviour of the functions Ck{x) 
From the first equation it is seen that 



Co(x) = -[2mV(x)]'/^ = -mujx(l + V /.a;^/^ = ^y C^x' (11) 

where the minus sign is chosen from boundary conditions, and coefficients Cf are defined 
by parameters of the potential as 

C',=-mu; C^ = ^(£c'^Cl^-2mf,) t>l. (12) 

Because the point x = is a simple zero for the function Co{x) , the function Ck{x) has 
a pole of the order of {2k — 1) at this point and consequently can be represented by a 
Laurent series 

oo 
1=0 

Then, according to the theorem of residues, the quantization conditions (10), 
expressed in terms of the Laurent series coefficients, take an especially simple form 

C2k-2 = ^Kk (14) 

where the symbol 5i^k is the Kronecker delta. 

Substituting the expansions (11) and (13) into (6) and equating coefficients of equal 
powers of a;, we derive 

(3 - 2A; + {)Ct^ + E E Ci,Ct^ = -2mEk6,,2k-2. (15) 

j=Op=0 

Sorting out the case of i ^ 2k — 2 yields the recursion relation for obtaining the 
coefficients Cf : 

Cf = -^[(3 -2k + t)Ct' + E E Cl,C^^ + 2 E C;cl^] (16) 

^"-"O j=l p=0 p=l 

whereas putting i = 2k— 2 we would find the recursion formula for the energy eigenvalues 



2mE, = -C^,}, - E E qC'2\-i2„, . (17) 

j=0 p=0 

Thus, equations (16) and (17) determine coefficients of perturbation expansions of 
energy eigenvalues and eigenfunctions in the same simple form both for the ground and 
excited states. 
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It should be noted, that just conditions (9), represented in the form (14) simphfy 
the consideration of the perturbation theory by means of semiclassical expansions. It 
becomes more evident under comparison of our technique with the rescaled version of 
the WKB approach [17-19], where the order of the energy is taken into account by the 
equality E = he. In this case, the coefficients Ck{e,x) in equation (7) have poles, too. 
The WKB quantization condition (8) then reads as 

oo 

J2f^'"^ResCk{€,x)=n (18) 

fc=0 

with residues polynomial with respect to e. By truncating the series in the left-hand side 
we arrive at the equation for determination of approximations to the energy eigenvalues 
e, and only after subsequent re-expanding, we can restore the results of perturbation 
theory [17-19]. 

4. Discussion and examples 

From equation (17) it is readily seen that for the energy eigenvalues, when k = 1, we 
immediately have the oscillator approximation 

Ei = uj{n+^) (19) 

and with k = 2 one obtains the form familiar from standard textbooks [15] 

^--4!S("^ + " + i)-^J^("^ + " + ^'- P°' 

It is easy to demonstrate that in the case of the harmonic oscillator our technique 
restores the exact solution for the wave functions as well. 

Putting, for simplicity, h = m = u = 1, from equations (12), (13) and (16) we find 

Co(x) = -X Ck{x) = dkX^-^'' A; > (21) 

where 

fc-i 
2dk = (3 — 2k)dk_i + y^ djdk-j k > 1 

U (22) 

d\ = n. 

Carring out the integration of the function Cq{x) we obtain the exponential factor of 
the eigenfunction. Its remaining part is a polynomial that satisfies the equation 

oo 

PjPn = J2dkx'-'' (23) 

k=l 

and, consequently, has a form 

mo 

Pn{x) =x''Y. Oia;^' CT = or 1 n = 2mo + a. (24) 

j=0 



On the basis of equation (23), the polynomial coefficients , a^, are determined by the 

system 



(n - 2m - a)am + c^sOm+i + ••• + dmo-m+iann, = 0. (25) 



The combination of these equations, multiplyed by a suitable dj with a view of taking 
into account equation (22), arrives at the following relation between two consecutive 
coefficients: 



(2m + (T + 2)(2m + a+l) 
4(mo — rn) 



that is the recursion formula for the Hermite polynomials ( see, for instance [30]). 

And at last, as an example, we consider the anharmonic oscillator with the potential 



V{x) = \x^ + \\x^. (27) 



Though this oscillator is widely discussed in the scientific literature [31-33], the analytical 
expressions for the perturbation expansions of energy eigenvalues presented in [31] are 
incorrect. The correct expansion coefficients obtained by means of formulae (16) and 



(17) have the form 

E3 = 5 A 2-^{4:n^ + 6n^ + 8n + 3) 

E5 = -A^ 2"^ (1572 n^ + 3930 n^ + 12220 n^ 
+14400 n^ + 11528 n + 3495) 

E7 = 5 A^ 2"^^ (23592 n^ + 82572 n^ 

+418236 n^ + 839160 n^ + 1523968 n^ 
+1488078 n^ + 939884 n + 247935) 

Eg = -\^2'^^ (45804660 n^ + 206120970 n^ 



+1471569960 n^ + 4188597000 n^ (28) 

+12317818548 n^ + 20804002800 n^ 
+29394281120 n^ + 25244303400 n^ 
+13898196592 n + 3342323355) 
Eu = 5X^2-^^ (1023655464 n" + 5630105052 n^° 
+52379661180 n^ + 193482687420 n^ 
+801071289576 n^+ 1940241040920 n^ 
+4424265058200 n^ + 6633369121920 n^ 
+8108461519360 n^ + 6378900376878 n^ 
+3214574914460 n + 725076383025) . 

In conclusion, a new useful technique for deriving results of the logarithmic 
perturbation theory has been developed. Based upon the h -expansions and suitable 
quantization conditions, new handy recursion relations for solving the bound-state 
problem for an anharmonic oscillator within the framework of the one-dimensional 
Schrodinger equation have been obtained. Avoiding the disadvantages of the standard 
approach these formulae have the same simple form both for ground and exited states 
and provide, in principle, the calculation of the perturbation corrections up to an 
arbitrary order in the analytical or numerical form. The extention on the three- 
dimensional case and the relativistic equations will be published elsewhere. 
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